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Abstract. A linear system of differential equations describing a joint motion of a 
thermoelastic porous body with a sufficiently large Lame's constants (absolutely 
rigid body) and a thcrmofluid, occupying porous space, is considered. The rig- 
orous justification, under various conditions imposed on physical parameters, is 
fulfilled for homogenization procedures as the dimensionless size of the pores tends 
to zero, while the porous body is geometrically periodic. As the results we derive 
Darcy's system of filtration for thermofluid, depending on ratios between physical 
parameters. The proofs are based on Nguetseng's two-scale convergence method of 
homogenization in periodic structures. 
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Introduction 

In the present publication we consider a problem of a joint motion of thermoelas- 
tic deformable solid (thermoelastic skeleton), perforated by a system of channels 
(pores) and incompressible thcrmofluid occupying a porous space. In dimensionless 
variables (without primes) 

x' = Lx, t' = Tt, w' = Lw, 9' = 9, — poPs, Pf — Pop/, 

the differential equations of the problem in a domain e M"^ for the dimensionless 
displacement vector w of the continuum medium and the dimensionless temperature 
9, have a form: 

(1) °^^P~g^ ^ diV:rIP, 

99 d 

(2) "■'"^^ at ^ div^(aKVa;6') - a^— div^jU;, 



(3) P = xa^Ji[x, — j + (1 - x)aAB(x, w)^{q + ae,(l - x)d% 

(4) q = P+ — l^+(^ofX9, 

ap ot 

(5) p + xo^pdiY^w = 0. 
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The problem is endowed with initial and boundary conditions 

Buj 

(6) w\t^o=0, ^|t=o = 0, 0|t=o=O, xen 

at 

(7) w^wo, 0^00, xeS = dn, t > 0. 
Here and further we use notations 

P = XPf + {'^- x}Ps, Cp = xcpf + (1 - x)cps, 
Ok = x<^Kf + (1 - x)o:ks, ae = xaef + (1 - x)<^es- 
In this model the characteristic function of the porous space x(^) is a known func- 
tion. 

For more details about Eqs.([T|)- ([5]) and description of dimensionless constants 
(all these constants are positive ) see [1]. 
We accept the following constraints 

Assumption 0.1. domain ft = (0, 1)^ is a periodic repetition of an elementary cell 
— eY , where Y — (0, 1)^ and quantity 1/e is integer, so that il always contains 
an integer number of elementary cells Yf . Let Ys be a "solid part" of Y , and the 
"liquid part" Yf ~ is its open complement. We denote as 7 = dYf D dYg and 7 
is -surface. A porous space fl'y is the periodic repetition of the elementary cell 
eYf , and solid skeleton f2^ is the periodic repetition of the elementary cell eYg. A 
boundary — dfll D dQ'y is the periodic repetition in f2 of the boundary £7. The 
"solid skeleton" fig is a connected domain. 

In these assumptions 

Xix) ^X^x) ^xix/e), 
Cp = c^ix) = x^{x)cpf + (1 - x%x))cps, 
p = p<'{x) = x^{x)pf + (1 - X^(a;))Ps, 
= otK^x) = xHx)a^f + (1 - x^(a;))a^s, 
ae = agix) = xHx)aef + (1 - xHx))aes, 
Let £ be a characteristic size of pores I divided by the characteristic size L of the 
entire porous body: 

I 

Suppose that all dimensionless parameters depend on the small parameter e and 
there exist limits (finite or infinite) 

lima^(£)=^o, limaA(£) = Ao, lima^(£)=ro, 

eXO e\0 6\0 

lim a,f{e) — 770, hm ap{e) — p*, hm ai,{e) — vq, 

£\,0 £\0 £\,0 

lim a^f = Ko/, lim a„s(£) = kqs, lim aef{e) = Pof, 

s\0 e\0 e\Q 

lim ags{e) = Pos, lim — ^ = /ii. 

£\,0 £\0 £^ 

We restrict our consideration by the case when when 

/io = 0; Ao ?7o = 00; Pos, /?o/, J^o tq < 00, 
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< ^1 + To, p*, Kof, KOs < OO. 

The condition Aq = ryo = oo means that the sohd skeleton is an absolutely rigid 
body. 

Using Nguetseng's two-scale convergence method [2J [Sj we show that the lim- 
iting regime as e tends to zero is described by two types of Darcy's like system 
of equations of filtration for the velocity of the liquid component, coupled with 
corresponding heat equation. 

Different isothermic models have been considered in [4j-[T0j. The first research 
with the aim of finding limiting regimes in the case when the skeleton was assumed 
to be an absolutely rigid isothermal body was carried out by E. Sanchez-Palencia 
and L. Tartar. E. Sanchez-Palencia O Sec. 7.2] formally obtained Darcy's law 
of filtration using the method of two-scale asymptotic expansions, and L. Tartar 
[3 Appendix] mathematically rigorously justified the homogenization procedure. 
Using the same method of two-scale expansions J. Keller and R. Burridge [4] derived 
formally the system of Blot's equations in the case when the parameter was of 
order e^, and the rest of the coefficients were fixed independent of e. Under the 
same assumptions as in the article [4], the rigorous justification of Blot's model 
has been given by G. Nguetseng [6] and later by A. Mikelic, R. P. Gilbert, Th. 
Clopeaut, and J. L. Ferrin in [8l[9l[7]. The most general case has been studied in 

m- 

§1. Formulation of the main results. 

As usual, equations ([I])-(l2l) are understood in the sense of distributions. They 
involve the equations ([T])- (l5|) in a usual sense in the domains (7^ and and the 
boundary conditions 

(8) [d] = 0, [w] =0, xoe T', t > 0, 

(9) [P] = 0, [a%v,0] = 0, Xoe t > 

on the boundary , where 

[(p]{xo) = f{s){xo) ~ (p(^f){xo), 

'P{s){xo) ^ lim (p{x), Lp(j){xo) ^ lim ip{x). 

X — » XQ X — > XQ 

X e ci^ X e Q^j: 

There are various equivalent in the sense of distributions forms of representa- 
tion of equations ll])-© and boundary conditions (l8])-(l9l). In what follows, it is 
convenient to write them in the form of the integral equalities. 

Definition 1.1. Functions (w^ , 9'^ , , are called a generalized solution of the 
problem ([l}- ([7]), if they satisfy the regularity conditions in the domain VIt — 

n X (o,T) 

(10) w", B){x, u;^), div^w", g^ p^ 9^,V^0'' G L'^i^T) 

in the domain 57^ — ^ x (0,T), boundary conditions ^ with functions 

equations 

(11) q'=f + -% + aefX' 

Up ot 

(12) + x'^ctpdiv^w^ = 



ene 
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a.e. in ^t, integral identity 

(13) Jqt^ 

{(1 - x')axHx, w') - [q' + ag.il - x'W)l} ■ ©(x, t^)) dxdt = 

for all smooth vector-functions ip — ^p{x, t) such that <p|af2 = '^\t=T = / dt\t=T — 
and integral identity 

(14) j {{arcie^ + aldiv^w')^^-alV^B' -V^i + ^i^dxdt^Q 

for all smooth functions ^ = ^(a;, t) such that £\dn — C|t=T — 0. 

In (|13p by ^ : i? we denote the convolution (or, equivalently, the inner tensor 
product) of two second-rank tensors along the both indexes, i.e., A : B = {B* o 

^) ~ J2i,j=i AijBj^. 

We suppose the next assumption to be held: 

Assumption 1.1. Sequences {{,JaJ^+,JaZ)x^\^ dui'^^^/ dt\} , {,/a\{\~x^)^ dwQ/ dt} , 
{{^/a^+ ^/a^)divxdwQ/dt}, {y/a^d'^WQ/ dt^} , {^/a^dOQ/ dt} , {V^qI '^^^ uniformly 
in e hounded in L^{fl). 

In what follows all parameters may take all permitted values. If, for example. 
To = , then all terms in final equations containing this parameter disappear. 
The following Theorems 1 1 . ll - [TT2l are the main results of the paper. 



Theorem 1.1. For all e > on the arbitrary time interval [0,T] there exists a 
unique generalized solution of the problem ([T]) - ([7]) and 

(15) ^ina^J\w%t)\+^x'\^.w^t)\ + il-x')V^\y.w%t)\\\2,n<Co, 

(16) Wh,n^ + ||V,r|l2,a, < Co, 

(17) Uh.nT + \\p'\\2,n^ + — ll^l|2,o, < Co 

ap at 

where Co does not depend on the small parameter e. 

Theorem 1.2. Sequences {x'^w^}, {p^} and {q^} converge weakly in LP'{Q.t) to 
, p, and q respectively. The sequence {O"^} converge weakly in L^{{0,T); W2^(f2)) 
to function 6. Functions admit an extension from f2g x (0,T) into such 
that the sequence {u^} converge strongly in L^((0, T); W2^(0)) to zero. Weak and 
strong limits v = dw^ /dt, q, p and 6 satisfy in the domain fir the state equation 

(18) q = p + i'oP^^dp/dt + mPofO, 
the continuity equation 

(19) {l/p^)p+ div.xW ^0, 
the heat equation 

(20) r„c;^-^| = rfz..(i3«.V^) 

ot p* ot 
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/ Bi{ni,t - t) ■V^q{x,T)dT 



and Darcy 's law in the form 
(21) 

in the case of tq > and /ii > 0; Darcy's law in the form 

(22) -u = -B2(mi) • 

in the case of tq — and, finally, Darcy 's law in the form 

1 /■* 

(23) v = (ml-Sa)-/ V^q{x,T)dT 

Top/ Jo 

in the case of ni — 0. The problem is supplemented by the boundary condition 

(24) {v ~ vo) ■ n{x) = 0, x e S, t>0 

for the velocity of the liquid component and boundary and initial conditions 

(25) To6'(cc, 0) = 0, 9{x,0) = 9oix) 
for the temperature. 

In Eq. pop symmetric strictly positively defined matrix is given below by for- 
mula (|53|) and in Eqs. (|2i p - (|23|) n{x) is the unit normal vector to S at a point 
X (z S, and symmetric strictly positively defined matrices i?i(/ii, t), i?2(/^i), and B^ 
are defined in jlOj . 

§2. Preliminaries 

2.1. Two-scale convergence. Justification of Theorems I1.1H1.2I relies on 
systematic use of the method of two-scale convergence, which had been proposed 
by G. Nguetseng [2] and has been applied recently to a wide range of homogenization 
problems (see, for example, the survey 0]). 

Definition 2.1. A sequence {(p^} C L'^{nT) is said to be two-scale convergent to a 
limit ip e L^{flT X Y) if and only if for any 1-periodic in y function a = crix, t, y) 
the limiting relation 

(26) lim / (f'^{x,t)a {x,t,x/e) dxdt = / / ip{x,t,y)a{x,t,y)dydxdt 
holds. 

Existence and main properties of weakly convergent sequences are established 
by the following fundamental theorem [51 [3] : 

Theorem 2.1. (Nguetseng's theorem) 

1. Any bounded in L^{Q) sequence contains a subsequence, two-scale convergent 
to some limit ip S L'^(Q,t x Y). 

2. Let sequences {v?^} and {eS/ x^^} be uniformly bounded in L'^{VLt). Then there 
exist a 1-periodic in y function (/? = ip{x,t,y) and a subsequence {^p'^} such that 
ip^S/yip e L^{flT X Y), and ip^ and e'VxV'^ two-scale converge to (p and ^yP>, re- 
spectively. 

3. Let sequences {(p*^} and x'P^} be bounded in L^{Q). Then there exist func- 
tions (p> G L'^{^It) and G L^{ilT x Y) and a subsequence from {(p^} such that 

is 1-periodic in y, 'S/yip € L^iVlj' x Y), and pf and 'S/x'P'^ two-scale converge to tp 
and V x'pix, t) + V yil:{x, t, y), respectively. 
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Corollary 2.1. Let a e i^(i^) a-nd <t^{x) a{x/e). Assume that a sequence 
{v?'^} C i^(f2T) two-scale converges to if £ L^{Q,t x Y). Then the sequence u'^ip^ 
two- scale converges to crLp. 

2.2. An extension lemma. The typical difficulty in homogenization problems 
while passing to a limit in Model {NBY as e \ arises because of the fact that the 
bounds on the gradient of displacement \I may be distinct in liquid and rigid 
phases. The classical approach in overcoming this difhculty consists of constructing 
of extension to the whole Vl of the displacement field defined merely on fi^. The 
following lemma is valid due to the well-known results from [!![ I12j . We formulate 
it in appropriate for us form: 

Lemma 2.1. Suppose that assumptions of Sec. 1.2 on geometry oj periodic struc- 
ture hold, e W2{i^s) '^''^'^ ip^ = on Sg = dfl^ H dft in the trace sense. Then 
there exists a function cr^ € W2 {^) such that its restriction on the sub-domain fl^ 
coincide with ip^ , i.e., 

(27) (i--x'{x)W(x)-r{x)) = o, xen, 

and, moreover, the estimate 

(28) Ik^|l2,n < C||7A^||2,o;, ||V,a^||2,o <q|V,^^||2,f2; 

hold true, where the constant C depends only on geometry Y and does not depend 
on e. 

2.3. Priedrichs Poincare's inequality in periodic structure. The fol- 
lowing lemma was proved by L. Tartar in [5, Appendix]. It specifies Friedrichs- 
Poincare's inequality for e-periodic structure. 

Lemma 2.2. Suppose that assumptions on the geometry offl^f hold true. Then for 

o 

any function (p &W2 (f^j) the inequality 



[ \ip\^dx < Ce^ [ \\7x^\^dx 



(29) 

holds true with some constant C, independent of e. 

2.4. Some notation. Further we denote 

1) = Iy ^dy, {^)yj = '^dy, = Iy, ^dy. 

2) If a and b are two vectors then the matrix a ® b is defined by the formula 

(a (g) b) • c = a(b • c) 

for any vector c. 

§3. Proof of Theorem HH] 

Under restriction tq > estimates (fT5|) - (fT6|) follow from estimates 



max (VoaIKI - x'')^ xW^ {t)\\2M + ^/a^\\^T-{t)\\2,n- 



^\\x'divxw'{t)\\2,n + V^\\0'{t)\\2,n) + \\VxO'\\2,nr + 

(30) y/aJ^Wx""^ x-Q^\\2,nT + \/a^llx^div^-^||2,nT < C'o, 

where Co is independent of e. These estimates we obtain if we multiply equation 
for w'^ by dw'^/dt, equation for 9'^ multiply by 0^, integrate by parts and sum the 
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result. The same estimates (|30p guaranties the existence and uniqueness of the 
generahzed solution for the problem (IT])-(l7]). 

Estimates ((T7| for pressures follows from Eas. (fTT1) - P^ and estimates (PU)) . 

Estimation of w'^ and 9^ in the case tq = is not simple, and we outline it in 
more detail. 

First of all we use estimates (l30l) in the form 



(31) V^\\x''^.w'{t)h.n + V^IKl - x')^xw^t)h.n < Co- 

Next, on the strength of Lemma 12.11 we construct extensions and ttg of the 
functions and Wq from fl^ into fi^-, such that u'^ — and Uq = Wq in fi^, Uq 
uniformly with respect to e bounded in W2 (f^) and 

\\u' - uU2.n < C\\V,{u' - u^h.n < ^||(1 - x")\/^V,(-u;" - w'o)hM- 



After that we estimate ||ii'^||2,n with the help of Friedrichs-Poincare's inequality 
in periodic structure (lemma [2?2)) for the difference {u^ — w^) and estimates (|3T|) : 

Ww'h.n < Wu'hs^ + \\u' - w'\\2,n < \\u'h,n + Ce\\x'V^{u' - w')\\2,n 
< ||m1|2,o + Ce||V^M^||2,o + Ciea;,^)\\x'y/^Vccw'\\2,n < 



Ce -1 
|m^||2,o + ^=||\/aIV^t(;''||2,o| + C{eaf,^)\\y^V^w''\\2,m < Co- 
Vox " ' 



The norm ||0'^||2,f2 we estimate with the help of the usual Poincare's inequality 
for the difference {6^ — (^q) and estimate (|30|) . 

The rest of the proof is the same as for the case tq > 0. 

§4. Proof of Theorem HH 

4.1. Weak and two-scale limits of sequences of displacement, tem- 
peratures and pressures. On the strength of Theorem ll.il the sequences {9^}, 
{V0^}, {p^}, {q^}, and {w^} are uniformly in e bounded in L'^{^It)- Hence there 
exist a subsequence of small parameters {e > 0} and functions 6, p, 5, tt and w 
such that 

r 9, Vr ^ Ve, p' -> p, q" ^ q, 

weakly in L^(JIt) as e \ 0. 

Due to Lemma [^?T1 thcrc is a function u'^ E L°°((0, T); W2{fl)) such that u.^ — w'^ 
in ns X (0, T), the family {u"} is uniformly in e bounded in L°°((0, T);W:^{n)) and 

u^^O strongly in L2((o,T);VF2^(r2)) 

as e \ 0. 
Moreover, 

(32) x'cx^.Hx,w')^0, 

strongly in L^{nT) as e \ 0. 

Relabeling if necessary, we assume that the sequences converge themselves. 

On the strength of Nguetseng's theorem, there exist 1-periodic in y functions 
Q{x,t,y), P{x,t,y), Q{x,t,y), and W{x,t,y), such that the sequences {\79^}, 
{p"^}; and {w^} two-scale converge to \7x9+\/yQ{x, t, y), P{x, t, y), Q{x, t, y), 
and W{x,t,y) respectively. 

Note that the sequence {div^ju;'^} weakly converges to divxW. 
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The same arguments we apply for the functions Wq and ^q- 

Wq Wo 

weakly in L'^{flT) as e \ 0, 

&o ^ ^0 

weakly in i^((0, T); W2{^)) as e \ and the sequence {(Hv^Wq} weakly converges 
to divxWQ. 

4.2. Micro- and macroscopic equations. 

Lemma 4.1. For all x e ^ and y (z Y weak and two-scale limits of the sequences 
{^'^}; {tt^}; {q^}, {W^}, {V:ri^^} and {S/xW^} satisfy the relations 

(33) P = (l/m)xp, Q = il/m)xq = (l/m)x(p + vop-^dp/Ot) + PofxO; 

(34) (l-x)W^ = 0; 

(35) q = P + vop'^dpl dt + mf3ofO; 

(36) (l/p,)p + diVxW = 0; 

(37) ({w - Wo) ■ n{x) = 0, x e S; 

(38) diVyW = 0. 

Proof. The weak and two-scale limiting passage in Eq.([TT|) yield that Eg. ([551 and 

(39) Q{x,t,y)^P{x,t,y) + ^dPix,t.,y)/dt + riy)9{x,t), yeYf. 

p* 

In order to prove equation ((33|) . into Eq. (fT3|) insert a test function ■0^ = ei/j {x, t,x/e), 
where xp{x,t,y) is an arbitrary 1-periodic and finite on Yf function in y. Passing 
to the limit as £ \ 0, we get 

(40) VyQ - 0, ye Yf. 

Next, fulfilling the two-scale limiting passage in the equality 

(1 - X')P' = 

we arrive at 

{1-X)P = 

which along with Eas. (P5)l - (HDl) justifies equation ([55]) . 

Eq. (|36p and boundary condition p7p are derived quite similarly if we represent 
Ea. (fT2|) in the form 



(41) — -I- diva;(u;^ - wD = (1 - x^)diva:U'' - div^wL 

dp 

multiply by an arbitrary function ■0'^ = ip{x,t), integrate, and then pass to the 
limit as e \ 0. Using now in ((4T|) test functions in the form xp^ = eif>{x, t,x/e) we 
obtain (|38l) . 

In order to prove Eq. ([M]) it is sufficient to consider the two-scale limiting relation 

in 

□ 
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Lemma 4.2. For all {x,t) G JIt the microscopic equation 
(42) diVy{K{y){V.J + VyQ)} = 0, y G F 

holds true. Here K = KQfX + kos(1 — x)- 



Proof. First of all, using continuity equation p2|) we rewrite the heat equation in 
the form 

(43) arc;^ = dW^ialV^e^) - a«.(l - x^)|(div.«-) + 

Substituting now a test function of the form ip'^ — etp {x,t,x/e), where xlj{x,t,y) 
is an arbitrary 1-periodic in y function vanishing on the boundary 9ri, into cor- 
responding integral identity, and passing to the limit as £ \ 0, we arrive at the 
desired microscopic relation on the cell Y . □ 

In the same way, using a test function independent of the fast variable y/e, we 
get from Eq. (|i5)l 

Lemma 4.3. For all {x,t) G fir the macroscopic equation 

(44) roCp|^ - = dtv.iikoV,^ + (ifVe)y} 

hold true. Here Cp = mcpf + (1 — m)cps, kq = {K)y- 

Now we pass to the microscopic equations for the velocities in the liquid. 
Lemma 4.4. Let V = xdW/dt. Then 
dV 

(45) ToPf-^ ^ fiiAyV -^yR^^^q, yeYf, t > 0; 

(46) V = 0, yG7; V(y,0) = 0, y G F/ 

in the case fii > 0, and 

dV 

(47) ToPf-g^ ^-VyR-y,q, yeVf, t>0; 

(48) y-n = 0, yG7; V(y,0) = 0, y G F/ 

in the case /ii = 0. 

In j?Q. (|48|) n is the unit normal to 7. 

Proof. Differential equations ((45|) and ((47|) follow as e \ from integral equality 
(fT3| with the test function ■0 = cp{xe^^) ■ h{x,t), where tp is solenoidal and finite 
in 17. 

Boundary condition in ([46]) is the consequence of the two-scale convergence of 

a_ 1 

{a^Wxw'^} to the function plWyW{x,t,y). On the strength of this convergence, 
the function \/yW{x,t,y) is L^-integrable in Y. The boundary condition (pS)) 
follows from Eq. ((38| . □ 



4.3. Homogenized equations. 
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Lemma 4.5. Weak and strong limits q and 9 satisfy the initial-boundary value 
problem 

(49) roCp^-^^ = div4B''-V^), x e n, t > 0; 

at ot 

(50) To9ix, 0) = 0, xen; 9{x, 0) = 9oix), xeS,t> 0, 

where a symmetric strictly positively defined matrix is defined by formula (j53p . 
Proof. For i — 1,2,3 we consider the periodic in y model problems 

(51) diYy{K{y)[VyQ, + e,)} = 0, y G F, 

where (ei, 62, 63) are the standard Cartesian basis vectors, and put 

3 

(52) e = (E> e,) • V^,9. 

1=1 

Then Q solves the problem (|5T|) and 

3 

(53) B*^ = kol + Y^iK'^y^Dr <E)e,). 

All properties of the matrix B^ arc well known ( see [S], [12] ). □ 

Lemma 4.6. The strong and weak limits 9, v = {V)y; , P and q satisfy in il-r 
equation (j36p , Darcy 's law in the form 



(54) 



V = - Bi{fii,t-T) ■ \7xq{x, T)dT 
Jo 



in the case of Tq > and fii > 0, Darcy 's law in the form 

(55) v = ~B2ifii)-V^q 

in the case 0/ tq — and, finally, Darcy 's law in the form 

(56) (ml-Sa)-/ V^q{x,T)dT 

ToPf Jq 

in the case of fii = 0. The problem is supplemented by the boundary condition 

(57) {v - vo) ■ n{x) = 0, x e S, t > 0, 

In Eqs. (|54p - (j57p n{x) is the unit normal vector to S at a point x Cz S , and 
symmetric strictly positively defined matrices Bi{pi,t), B2{fJ-i), and B3 are given 
in [TU]. 

The proof of these statements repeats the proof of Lemma 5.8 in [TU] , 
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